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ABSTRACT 

The cloud-in-cloud problem is studied in the context of the extension to non-Gaussian 
density fields of the Press-Schechter approach for the calculation of the mass func- 
tion. As an example of a non-Gaussian probability distribution functions (PDFs) we 
consider the Chi-square, with various degrees of freedom. We generate density fields 
in cubic boxes with periodic boundary conditions and then determine the number of 
points considered collapsed at each scale through an hierarchy of smoothing windows. 
We find that the mass function we obtain differs from that predicted using the Ex- 
tended Press-Schechter formalism, particularly for low values of a and for those PDFs 
most distinct from a Gaussian. 
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1 INTRODUCTION 

One of the distinctive features in the Universe is the presence 
of gravitationally collapsed structures, like galaxies, groups 
and clusters of galaxies. The distribution of masses of these 
structures, usually called the mass or multiplicity function, 



has been determined observational ly (|Ashman, Salucci k 



Persic 1993| ; [Henry fe Arnaud 199l| ; [Eke ct al. 1998f ; jMarke-, 
vrtch"l*998[ ) , and is one of the most important characteristics 



of the Universe that proposed cosmological models attempt- 
ing to explain the formation of structure need to reproduce. 



(but see Lucchin & Matarrese 1988) applied in the context 
of structure formation models where the perturbations in- 
duced in the density field have a Gaussian random-phase 
distribution independently of the scale considered. This as- 
sumption is not only the simplest to take, but is also ex- 
pected when perturbations are produced by an inflation- 
ary phase in the very early Universe (e.g. Liddle & Lyth 
1993). However, there is at present a renewed interest in 
structure formation models which predict a non-Gaussian 
density distribution, either within the context of inflation 
(IPeebles 19991 ; [Salopek 1999] |Martin, Riazuelo fe Sakellar- 



In order to proceed with a comparison between these ob- 
servations and the theoretical expectations of different strnc- 



ture for -nation models, it is of fnndn.nwitnl inmortn.nen tn he 



199E 



iadou 199S|), or as a result of the dynamics of topological 
defects (| Avelino et al. 1998a ; Albrccht, Battye fe Robinson 



Avelino, Caldwell & Martins 1999 



able to lpredict with reasonable accuracy the mass function 



Contaldi, Hind- 



associated with the theoretical models. The most direct way 
of doing this is to perform N-body simulations (for a recent 
review see Bertschinger 1998), where a distribution of dark 
matter particles alone, or in conjunction with gas particles, 
is evolved under gravity. However, these simulations take 
a considerable time to complete, being impracticable if a 
large number of structure formation models are being stud- 
ied simultaneously. Another method for estimating the the- 
oretical mass function is to use analytical approximations. 
Among the several that have been proposed, the framework 



marsh & Magueijo 1999). In an attempt to compare their 
predictions with the observed mass function, particularly 
at the scale of galaxy clusters (Chiu, Ostriker & Strauss 
1998; Koyama, Soda & Taruya 1999; Robinson, Gawiser & 
Silk 1999a,b; Willick 1999) a particular generalization of the 
Press-Schechter framework has been used, so that the as- 
sumed density field no longer needs to be Gaussian. This 
so-called Extended Press-Schechter (EPS) approach has re- 
cently been proved to be quite successful in reproducing the 
results for the mass function obtained from N-body sim- 



ulati ons with non-Gaussian conditions (Robinson & Baker 



put for ward by Press and Schechter (1974, hereafter PS) has 
proved the most successful in reproducing the mass function 
obtained through N-body simulations, albeit the very sim- 
plified assumptions that are made. 

Until recently the approach proposed by Press and 
Schechter to estimate the mass function was almost always 



199g ). 

The EPS method was obtained by closely following the 
reasoning behind the original PS work, in the hope that it 
would end up as successful in predicting the mass function. 
The fact that this seems to be the case increases our per- 
plexity as to why the general PS framework works at all, 
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given all the simplifications it entails, like spherical collapse 
and that at a given smoothing scale all the structures that 
form have equal mass. The particular issue we will study 
here is the so-called cloud-in-cloud problem, with the others 
addressed in forthcoming papers. This aspect of PS-based 
approaches to the calculation of the mass function has been 
investigated previously (Epstein 1983, 1984; Schaeffer & Silk 
1988; Peacock & Heavens 1990; Bond et al. 1991; Jedamzik 
1996; Yano, Nagashima & Gouda 1996; Monaco 1998), al- 
though always in the context of Gaussian initial conditions. 

Until now the cloud-in-cloud problem within the EPS 
framework has been dealt with in the same simplified man- 
ner as in the original PS derivation. In order to determine 
to what extent such a treatment of the problem is justi- 
fied, we follow the numerical approach laid down by Monaco 
(1997a,b; 1998), simulating density fields with non-Gaussian 
one-point probability distributions in cubic boxes with pe- 
riodic boundary conditions. 

In the first section we describe succinctly the Press- 
Schechter approach to the derivation of the mass function, 
and how it can be extended to accommodate non-Gaussian 
initial conditions. In the following section we discuss the 
cloud-in-cloud problem and how several authors have tried 
to avoid it, describing in detail the numerical approach we 
have chosen to study it. Finally, in the last two sections we 
present the results of our analysis and discuss their impor- 
tance to the proposed Extended Press-Schechter framework. 



2 EXTENDED PRESS-SCHECHTER 



The Press-Schechter theory was originally proposed (Pres: 



& Schetnter iyY4), m tne context ot initial Uaussian density 



perturbation^, as a simple analytical tool for predicting the 
mass fraction associated with collapsed objects with mass 
larger than some given mass thereshold M. This is obtained 
by measuring the fraction of space in which the evolved lin- 
ear density field exceeds a given density contrast 8 C , 



F(> M) 



O m (> M(R)) 



V(S)dS. 



(1) 



For spherical collapse in an Einstein-de Sitter universe <5 C 
equals 1.7, being almost insensitive to a change in the as- 
sumed background cosmology (e.g. Eke, Cole & Frenk 1996). 
In the case of Gaussian initial conditions then, 



F(> M) = ierfc ( 

1 ' 2 \V2a(R) 



(2) 



where erfc is the complementary error function and a(R) 
is the dispersion of the density field at the scale R. For a 
top-hat window, M is related to R via M = 47r7? 3 pb/3, 
with pb being the background density. The right hand side 
in expression (|^) is usually multiplied by a factor of 2, as 
originally suggested by Press and Schechter, so that 



dF(> M) 



(3) 



thus taking into account the accretion of material initially 
present in regions underdense at the smoothing scale R. The 
mass function of collapsed objects can be obtained simply 
by deriving expression (|^) (multiplied by the factor 2) with 
relation to M, and then dividing by pi,/M. 



The original PS approach can be easily generalized for 
non- Gaussian density perturbations, simply by considering 
in expression (Q) a non-Gaussian one-point probability dis- 
tribution function (henceforth PDF) V(S). In order for all 
the mass in the Universe to be accounted for the expression 
then needs to be multiplied by / = 1/ J" °° V{5)d5. This re- 
normalization is equivalent to multiplying by the factor 2 
in the Gaussian case. Surprisingly, this very simple exten- 
sion of the PS approach to non-Gaussian initial conditions 
was only taken seriously for the first time by Chiu, Ostriker 
and Strauss (1997). Since then it has been tested with some 
success for a few non-Gaussian structure formation models 



against N-body simulations (Robinson & Baker 1999) 



It should be noted that as long as the dispersion of the 
density field depends on the smoothing scale R considered, 
cr(R), the PDF for the matter density is necessarily scale- 
dependent, V(R,S). What is usually meant by saying that 
such PDF is invariant with scale, is that the reduced distri- 
bution, Vr(5) = V[R, S/a(R)]/a(R), is always the same, i.e. 
the shape of the reduced PDF does not depend on the scale 
under consideration. For example, the simplest inflationary 
models predict a Gaussian PDF for the matter density at 
all scales (e.g. Liddle & Lyth 1993), in the sense that at 
any scale the reduced form of the PDF is always equal to a 
Gaussian with zero mean and dispersion unity. 

In the case of alternative structure formation models 
one could have two other scenarios. In the more general 
one, the shape of the reduced PDF for the matter density 
depends on the smoothing scale being considered. This is in- 
deed what is expected when the matter density distribution 
is generated, for example, through the dynamics of topolog- 



ical d efects (Avelino et al. 1998b; Avelino, Wu & Shellard 



199S). The second possibility is that the reduced PDF for 
the matter density is scale-independent, but non-Gaussian, 
i.e. not equal to a Gaussian with zero mean and dispersion 
unity. Here, for simplicity, we will focus our study on this 
second scenario, with the knowledge that it can be easily 
generalized to the first one. 



3 SOLVING THE CLOUD-IN-CLOUD 
PROBLEM 

As we have previously shown, there is a fundamental dif- 
ficulty with the normalization of expression (Q). We have 
seen that in the Gaussian case a factor of 2 had to be intro- 
duced to correctly take into account the fact that material 
in initially underdense regions at some smoothing scale R is 
eventually accreted and incorporated into the collapsed ob- 
jects that form from the initially overdense regions at that 
scale. In other words, the factor 2 accounts for the mate- 
rial that although not in the regions predicted to collapse 
at the smoothing scale R, will nevertheless become part of 
collapsed objects associated with scales larger than R. Only 
by smoothing the density field at these scales would this ma- 
terial count as collapsed, but clearly this should happen in 
the first place when F(> M(R)) is calculated. This is the 
so-called cloud-in-cloud problem in the PS approach. 

Several authors have tried to find a more satisfactory 
solution for the cloud-in-cloud problem than just multiply- 
ing expression (Q) by 2, as proposed by Press and Schechter. 
The first to approach the problem were Epstein (1983, 1984) 
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Figure 1. For the reduced Gaussian PDF it is shown on the left F(a), and on the right its derivative with respect to cr. The solid lines 
correspond to the EPS prediction, while the dashed and dotted lines to the numerical results, respectively in the case of sharp-k and 
gaussian smoothing. 



and Schaeffer & Silk (1988). But it was only with the work 
of Peacock and Heavens (1990) and Bond et al. (1991) that 
a comprehensive framework was put in place to study the 
cloud-in-cloud problem and its ramifications. The assump- 
tions behind the two approaches are very similar, with the 
later using a more formal line of reasoning based on the the- 
ory of excursion sets. The method we will use here to study 
the cloud-in-cloud problem within the EPS framework was 
first used by Monaco (1997a,b; 1998) for the case of Gaussian 
density fields. It is basically a numerical implementation of 
the approach pioneered by Peacock and Heavens (1990) and 
Bond et al. (1991). Density fields with an assumed PDF, 
and which only differ in the scale at which the smoothing 
is applied, are generated in cubic grids. Starting from the 
largest scale, trajectories for all the points in the density 
field are then constructed with the values for the density 
contrast recorded at each smoothing scale. In the Gaussian 
case, and if the smoothing is performed with a sharp-k win- 
dow, the trajectories followed by the points are Brownian 
random walks. In this very particular example, the number 
of points which exceed some density contrast at any smooth- 
ing scale is equal to the number of points which though at 
such scale do not exceed the assumed density contrast, at 
some other larger scale do. In this case the re-normalization 
factor one needs to multiply expression (Q) with is exactly 
2, as Press and Schechter initially proposed. This was for- 
mally proved by Bond et al. (1991) using excursion set the- 
ory. Unfortunately, it is the only instance when their ap- 
proach can be used to solve the cloud-in-cloud problem. For 
any other smoothing window or probability distribution the 
random walk characteristic disappears, and the problem be- 
comes analytically intractable. In these cases, either one uses 
the method proposed by Peacock and Heavens (1990), and 
expanded in Monaco (1997a,b; 1998), which also has its lim- 
itations, or the numerical approach considered here. 



4 RESULTS 

Our numerical realizations of density fields were performed 
in cubic grids with 64 3 points. They were generated in the 
same manner as the density fields used to set initial condi- 
tions in N-body simulations (see e.g. Bertschinger & Gelb 
1991; Klypin & Holtzman 1997). We considered one exam- 



ple of a scale-independent reduced non-Gaussian PDF for 
the density contrast, the Chi-square, the PDF having been 
shifted so that its mean is zero. Such a PDF is expected 
in certain models of structure formation, involving either 
isocurvature density perturbations generated during an in - 
flationary period in the very early Un iverse ([Peebles 199Sh 



or cosmic string seeded pertu rbations ( Avelino et al. 1998b 



Avelino, Wu & Shellard 1999). We also considered the Gaus- 
sian case for comparison. 

The Chi-square has the added attractive feature that 
the reduced PDF remains approximately the same when 
smoothed on a variety of length scales using either gaus- 
sian (GAU) or top-hat (TH) windows (e.g. White 1999). 
The reduced PDF that results from the smoothing becomes 
increasingly different from the original as the number of de- 
grees of freedom, v, gets smaller, with the most important 
departures relatively to the original reduced PDF shape be- 
ing a decrease in the probability of S/a taking values around 
zero, i.e. the mean, and the appearance of a non-zero prob- 
ability of 5/cr taking values just outside the cut-off, equal 
to —u/y/2i/. Fortunately, these differences only become no- 
ticeable for v smaller than about 10. However, when the 
smoothing is performed using a sharp-k (SK) window, the 
departure of the smoothed PDF from the original increases 
considerably. Now, as soon as v becomes less than about 100 
the difference starts to show. Nevertheless, we will be only 
interested in the fraction of the reduced PDF that lies above 
5 c /a, and this part of the distribution is little changed by 
smoothing, even in the case of a SK window. 

With the above in mind, for different v values we as- 
sumed that the power spectrum associated with the den- 
sity contrast was a power-law, Vs{k) oc k n , where k is the 
wavenumber and n = —2, —1,0. We generated several sets 
of density fields for each combination of v and n values, 
such that all realizations within each set only differed in 
the smoothing scale applied. The three smoothing windows 
mentioned above were considered: sharp-k, top-hat, and a 
gaussian. However, we will only show the results for SK and 
GAU smoothing, given that those for TH smoothing turn 
out very similar to the GAU ones. Also, we opted for show- 
ing just the results obtained for n — —2, which is closest 
to the slope of t he matter power spectrum on the scale of 
galaxy clusters (Vlarkevitch 199c), as they are basically in- 
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Figure 2. The same as in Figure 1, but for the three reduced Chi-squarc PDFs considered in the paper. 



distinguishable from the results for the two other spectral 
indexes, n = — 1 and n = 0. 

On the left panel of Fig. 1 we show, for the reduced 
Gaussian PDF, the fraction of mass collapsed above a cer- 
tain smoothing scale as a function of the value of the disper- 
sion of the density field at that scale, F(a), calculated using 
the method presented here and the EPS approach. The same 
quantities are shown on the left panel of Fig. 2 for a reduced 
Chi-square PDF with v = 1, v = 10, v = 100 (note that the 
v = oo case is equivalent to a reduced Gaussian PDF). 

On the right panels of Fig. 1 and Fig. 2 we show 
dF(a)/da instead. The theoretical mass function can be eas- 
ily obtained by multiplying dF(a)/da by da/dM and divid- 
ing by pb/M. In the case of power-law spectra as the ones 
we consider here, er(M) oc M~ (n+3)/6 . 

In all cases the thereshold for the density contrast above 
which we assume a field point to be collapsed is S c = 1.7, 
which is equivalent to assuming spherical collapse in an 
Einstein-de Sitter cosmology. The results presented can be 
easily generalized for any other thereshold, S' c , by making 
the identification F'(a) = F(a * 5 C /S' C ). 

In the case of SK smoothing and in the limit of a re- 
duced Gaussian PDF, the EPS approach (which then simply 
reduces to PS) correctly predicts F(a), and thus dF(a)/da, 
as expected. The famous factor of 2 in the normalization is 
recovered. As the number of degrees of freedom decreases the 
assumed reduced PDF starts to differ from a reduced Gaus- 
sian, and the EPS prediction increasingly overestimates the 



numerically determined collapsed mass fraction and mass 
function. In the case of the GAU window (as mentioned 
before the results for TH smoothing are very similar), the 
numerical results deviate from those predicted through the 
EPS approach even when assuming a Gaussian PDF. For 
this particular PDF, we find that for small values of a (i.e. 
large mass scales), both the predicted collapsed mass frac- 
tion and the mass function seem to approach the result one 
would obtain if the PS approach was used without the nor- 
malization factor of 2. The same conclusion had already been 
reached analytically by Peacock and Heavens (1990), be- 
ing numerically confirmed by Bond et al. (1991) (see also 
Monaco 1997b). 



5 CONCLUSIONS 

We have shown that the EPS approach does not correctly 
take into account the existence of regions, which though at 
some smoothing scale are unable to pass a certain collapse 
thereshold are nevertheless able to do it at larger scales, 
in the estimation of the mass function. This problem was 
already known to exist within the PS framework, except 
when the smoothing of the density field was performed using 
a sharp-k window. Now, we have found that even when this 
window is used, the EPS approach cannot adequately solve 
the cloud-in-cloud problem. 

The mass function predicted by the EPS approach de- 
viates most from the numerical results for small values of 
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<t, which correspond to the large mass end in models where 
structure builds up hierarchically. The only N-body simu- 
lations that have been used to check whether the EPS ap- 
proach provides a good fit to the mass function were limited 
in range to values of a larger than around 0.5, except for 
a couple of points with relatively large error bars, and a t 
the limit of statistical significance (Robinson & Baker 1999). 



For these a values, the EPS prediction for the mass function 
has not yet entered the regime where it differs significantly 
from the numerical results obtained in this paper, particu- 
larly when one considers sharp-k smoothing. It is therefore 
not surprising that the N-body results seem to vindicate 
the EPS approach. However, our analysis should throw a 
note of caution. In the regime where the mass function is 
defined essentially through the abundance of rare density 



peaks, the EPS approach may not be reliable when one is 
dealing 1 with strongly non-Gaussian PDFs. This may affect 



some relcent conclusions regarding the gaussianity ot the den- 
sity field on large scales, drawn from the abundance of rich 
galaxy clusters at different redshifts and their present-day 
correlation length (Chiu, Ostriker & Strauss 1998; Koyama, 
Soda & Taruya 1999; Robinson, Gawiser & Silk 1999a,b; 
Willick 1999). It would be very interesting if larger N-body 
simulations could be carried out, able to extend the mass 
function further into the low-cr, rare events regime. 

The mass functions determined in this paper can still 
be further improved. Here we focused our attention on the 
cloud-in-cloud problem. Other issues were left untouched, 
more importantly, possible deviations from spherical col- 
lapse and the relation between smoothing radius versus the 
mass of the structures identified after each smoothing. We 
are presently looking at these issues. 
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